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Abstract 

We study the chromatic polynomials for m x n square-lattice strips, of 
width 9 < m < 13 (with periodic boundary conditions) and arbitrary length 
n (with free boundary conditions). We have used a transfer matrix approach 
that allowed us also to extract the limiting curves when n — > oo. In this limit 
we have also obtained the isolated limiting points for these square-lattice strips 
and checked some conjectures related to the Beraha numbers. 
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1 Introduction 



The antiferromagnetic g-state Potts model fl], 0, has many interesting features. 
First, for large enough q the model exhibits a nonzero ground-state entropy (with- 
out frustration). This is important as it provides an exception to the third law of 
thermodynamics @, |^. Second, the antiferromagnetic Potts model has a rich phase 
diagram that depends explicitly on the lattice structure (in contrast to the universal- 
ity typically enjoyed by ferromagnets). It therefore becomes interesting to elucidate 
which features of the lattice structure (e.g., coordination number, dimensionality, etc) 
correlate with the critical properties of the model (or the absence of criticality) . Fur- 
thermore, for each lattice G there is a number qc{G) such that for all q > Qc the Potts 
model is disordered at any temperature (including T = 0); at g = qdG) the system 
is disordered at all positive temperatures and has a zero-temperature critical point. 
Third, the study of Potts antiferromagnets is closely related to graph theory. Namely, 
the zero-temperature limit of the antiferromagnetic Potts model partition function is 
just the chromatic polynomial Pg{<i)'- 

limZG(g;T) = PG(g). (1.1) 

This quantity equals the number of ways of coloring the vertices of the graph G 
with q colors, with the constraint that no adjacent vertices have the same color 
ip. Finally, there are condensed-matter systems that can be modeled using Potts 
antiferromagnets. In particular, the insulator SrCr8_aGa4_2^^0i9 can be described 
using a three-state Kagome-lattice Potts antiferromagnet ^ ^. 

One interesting question is the value of the number qc{G) for the most common 
regular lattices G. There are general analytic bounds in terms of the coordination 
number A of the lattice (e.g., qc < 2 A ^), but they are not very sharp. An alternative 
approach, inspired by the Lee- Yang picture of phase transitions []TU[, is to study the 



zeros of the chromatic polynomial when the parameter q is allowed to take complex 
values (see Ref. |]ll| for a complete list of references). The best results concern 



families Gn of graphs for which the chromatic polynomial can be expressed via a 
transfer matrix of fixed size M x M: 

PgM = tr[A(g) r(g)"] (1.2a) 

M 
k=l 

where the transfer matrix T{q) and the boundary- condition matrix A{q) are polyno- 
mials in q, so that the eigenvalues {Xk} of T and the amplitudes {ak} are algebraic 
functions of q. 

The two best understood cases of the Potts antiferromagnet are the square and 
triangular lattices. In this paper we consider the square lattice. For q = 2 this model 
undergoes a finite-temperature second-order phase transition. For q = 3 there is 



strong analytic and numerical support for a zero-temperature critical point [|1J, |13 



HI , |T5|, [16], |T7|, |1^, 0, though a rigorous proof is (to our knowledge) still lacking. 
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For g = 4 Monte Carlo simulations show that the model is disordered even at zero 



temperature [|T9l. It thus seems clear that gc(sq) = 3. 

In the transfer matrix approach, the goal is to compute the transfer matrix T 
for a square-lattice strip of width m with some boundary conditions. With this 
matrix one can easily compute the chromatic-polynomial zeros of any finite strip 
m X n. In addition, one can also compute the accumulation points when n oo. 
According to the Beraha-Kahane- Weiss theorem |2^, ^ , the zeros when n —>■ oo can 
either be isolated limiting points (when the amplitude associated to the dominant 
eigenvalue vanishes) or form a limiting curve B (when two dominant eigenvalues 
cross in modulus). By studying the limiting curves for different values of the width 
m we hope to learn new features of the thermodynamic limit m ^ oo. 

The first results for the square lattice using this approach were obtained by Shrock 
and collaborators.^ They studied strips up to m = 6 (resp. m = 5) with cylindrical 
(resp. free) boundary conditions^ |]2^, They also computed square-lattice 



strips with cyclic and toroidal boundary conditions up to m = 4 p^ , ^ ^ p9 



50| . The results with free and cylindrical boundary conditions were extended to m < 8 



in Ref. |Tl|. The qualitative shape of the limiting curve found for these boundary 
conditions is quite similar to the one found by Baxter in case of the triangular lattice 
pl| , but with the zero-temperature critical point lying at 5'c(sq) = 3 rather than at 
gc(tri) = 4. 

However, there may be an important qualitative difference between the triangular- 
lattice and square-lattice cases, concerning the number of times that the infinite- 
volume curve B crosses the positive real g-axis. In physical terms this corresponds to 
the number of intervals of the real g-axis where the ground-state entropy 5*0 (g) has a 
distinct analytic form, or in other words to the number of "phases" that the model 



has as g is varied [^, 0. For many (but not all) finite strip widths m, we find that 
the limiting curve B crosses the real axis at a point go(m).0 It is reasonable to expect 
that there are infinitely many such widths m and that qo{iTi) tends to a limiting 
value go(oo) as m oo. One could be tempted to assume that go(oo) = qc] but 
this, it turns out, is not always true. Indeed, in the triangular-lattice case, Baxter's 



exact solution shows that go(tri) ^ 3.81967 while gc(tri) = 4; this is because the 
infinite- volume curve B crosses the real g-axis twice in the interval < g < oo. Thus, 
in general one cannot identify go with gc.0 On the other hand, our data for the square 

^ They actually used a generating-function approach which is equivalent to transfer matrices. 

^Let m (resp. n) denote the number of sites in the transverse (resp. longitudinal) direction of 
the strip, and let F (resp. P) denote free (resp. periodic) boundary conditions in a given direction. 
Then we use the terminology: free (mp x np), cylindrical (mp x np), cyclic (mp x np), and toroidal 
(mp X Tip). 

For those lattice strips where qq cannot be strictly defined because B does not cross the real 
axis, the limiting curve often includes a pair of complex-conjugate endpoints that are very close to 
the positive real q-axis. In these cases, we define Qq to be the endpoint closest to that axis with 
positive imaginary part. 

^ One can qualitatively think of Qc as the limit of the (complex) rightmost endpoints of the 
limiting curve B as the lattice strip width m increases. These endpoints can be defined in most cases 
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lattice suggest that qo{oo) may well be equal to qc = 3. If this is true, the qualitative 
shape of the infinite- volume curve B for the square lattice will be rather different from 
that found for the triangular lattice. In particular, region V2 in Figure 5 of ref. |^ 



will be "narrowed" along the g-axis so the points F (i.e., go) and C (i.e., Qc) coincide. 
Thus, region T>2 will be split into two complex-conjugate regions that may or may 
not intersect the positive real g-axis at qc- 

In this paper we extend the results of [|n[| by computing the transfer matrices 
for square-lattice strips of width 9 < m < 13 with periodic boundary conditions in 
the transverse direction and free boundary conditions in the longitudinal direction. 
We have chosen this type of boundary conditions because we expect it to provide 
a faster approach towards the thermodynamic limit than we would have obtained 
using free boundary conditions in both directions, while at the same time avoiding 
the technical complications implied by the imposition of doubly periodic (toroidal) 
boundary conditions The present computation has been possible thanks to a 



more efficient algorithm than the one employed in |TT|]; this improved algorithm allows 



us to handle transfer matrices as large as 498 x 498. The transfer matrices, of course, 
enable us to compute the chromatic polynomials Pmxn(g) for any finite length n. For 
9 < m < 11, however, we have been able to obtain the limiting curves B directly in 



the limit n 00, by using the same methods as in [|TT|. When possible, we have 
also computed the isolated limiting zeros for these square-lattice strips and checked 
whether the conjectures put forth in |Tl|, Section 7] hold true. These conjectures are 



related to the Beraha numbers |^ which seem to play a special role in the theory of 
chromatic polynomials (see Ref. [0] and references therein): 

S„ = 4cos2- = 2 + 2cos— , forn = 2,3,.... (1.3) 

n n 

We find that the limiting curves B for various system sizes are qualitatively quite 
similar; however, they show noticeable differences depending on whether the width 
m is even or odd. In particular, we conjecture that only the limiting curves for odd 
m cross the real axis, thus giving a well-defined value of go. Our best estimate for 
this quantity (assuming monotonicity with the strip width) reads 

go(sq) > 2.9161885031 (1.4) 

This estimate is quite close to the expected value for gc(sq) = 3. For the square-lattice 
strips with even width m we find for m > 8 that there is a small gap separating B 
from the real axis: the real part of the points closest to that axis is slightly smaller 
than the value of go for widths mil, and based on our numerical data we conjecture 
that the imaginary part of such points goes to zero as ^^'^ . 

We finally conjecture (based on our numerical evidence) that go(sq) = gc(sq) = 3. 
This equality has important implications on the analytic structure of the free energy 



for free boundary conditions in the longitudinal direction |T^, |3^, but not for periodic boundary 
conditions in the longitudinal direction p^ . When m grows large, those endpoints approach the 
real q-axis, and in the limit m = 00 they pinch that axis ai q = Qc- 
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in the thermodynamic hmit. In the zero-temperature triangular-lattice Potts antifer- 
romagnet, we have that goltri) 7^ Q'c(tri) and there are three domains of analyticity of 



the free energy ||3T[. Furthermore, all these three domains intersect the real g-axis; 
go and qc are the boundaries of such domains on the positive real g-axis. In the 
zero-temperature square-lattice Potts antiferromagnet, we have that go(sq) = gc(sq). 
Thus, on the real g-axis there are only two domains of definition whose boundary is 
precisely Qc- This is a valuable piece of information about the analytic structure of 
the free energy of this model whose solution is yet to be found. 

This paper is laid out as follows: In Section we explain the method used to 



compute the transfer matrices. The necessary background can be found in |jTT|. In 



Section ^ we expose our numerical results for square-lattice strips with cylindrical 
boundary conditions. Finally, Section ^ is dedicated to our conclusions. 



2 Transfer Matrix Algorithm 



The general theory of Potts model transfer matrices has been reviewed in |]TT 
Section 3], and in what follows we shall therefore limit ourselves to a rather concise and 
more informal description. The main emphasis is on the algorithmic improvements 
over Ref. |Tl|, which have enabled us to extend the computations to strips of width 



9 < m < 13 with cylindrical boundary conditions. 

Let G = {V, E) be a finite undirected graph with vertex set V and edge set E, and 
let {Je}e£E be a set of coupling constants. In the Fortuin-Kasteleyn representation 
[|37| , |38| the partition function of the g-state Potts model defined on G reads 

ZG{q,{ve})= E^'^""'^!!^- (2.1) 

E'CE e£E' 

where k{E') is the number of connected components (clusters) in the subgraph {V, E'), 
and Ve = exp(/3Je) — 1. The main advantage of this representation is that g can now 
be analytically continued to complex values. However, the price to be paid is that 
the Boltzmann weights contain the non-local factor q^^^'\ 

Although this non-locality would seem to inhibit the construction of the transfer 
matrix, the problem can be circumvented by working in a basis of connectivities Cm-, 
containing information about how the m spins in a given time slice t are interconnected 
in the portion of {V^E') that is prior to t. Clearly, the number Cm = \Cm\ of such 
connectivities is equal to the number of ways that m points on the rim of a disc can 
be connected in the interior of the disc, by means of a planar graph. The evaluation 
of these numbers is a standard combinatorial exercise and one arrives at the 
Catalan numbers, 

m\[m + Ij! 

A numerical implementation of the Potts model transfer matrix in the basis of 
Catalan connectivities was first described by Blote and Nightingale These au- 



thors also introduced a ranking of the connectivities, i.e. an injective mapping from 
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Cm to the set of integers {1,2,..., Cm}- Details on the construction of the inverse 
mapping were given in The idea is then to use the integers to index the entries 
of the transfer matrix, and the connectivities to work out the possible transitions 
between states along with their respective Boltzmann weights. When combined with 
a standard sparse matrix factorization, in which the entire transfer matrix is written 
as a product of elementary matrices each adding one edge of G, this algorithm is 
highly efficient, since each non-zero matrix element can be computed in time ~ m. 

Since we are considering specifically the zero-temperature antiferromagnet {ve = 
— 1 for all edges e E E) the dimension of the transfer matrix is actually much less 
than Cm-, since among the m points in a given time slice no pair of neighboring 
points can belong to the same cluster. It can be shown [0, Section 3.3] that the 
number Rm = \'R'm \ of antiferromagnetically allowed states TZm is given by the Riordan 
numbers, Rq = 1, Ri = Q and 

m-l [k/2\ 

i?^ = E(-ir"'"'E 27r^- (2.3) 
fc=0 j=0 ^ -^^ 

Furthermore, due to isotropy and the periodic boundary conditions, the connectivities 
TZm are invariant under the action of the dihedral group Dm on the cyclically ordered 
set of m points in the time slice considered. The appropriate basis states Sm are thus 
found from the Riordan states TZm by symmetrizing with respect to the elements 
of Dm- The resulting dimension of the transfer matrix SqCyl(m) = \Sm\ has been 
evaluated in |T1|, Table 2]. 



However, the sparse matrix decomposition does not respect any of these con- 
straints, and accordingly the number of connectivities in the intermediary states be- 
tween to subsequent time slices will still be Cm- To compute a given column T(g)|si) 
of the transfer matrix, with \s\) G Sm-, we therefore proceed as follows. First, each 
of the distinct states in the set {Dm\s\) : j = 1,2,..., 2m} C {Cm} is assigned a 
weight g'^d'^'^), where k{\si)) is the number of connected components in \si). Second, 
we act on this initial state with T(g), or rather with its sparse matrix decomposition. 
The final state obtained in this manner should again be symmetric with respect to 
all Dm , and it provides a non-trivial check of our algorithm that this is actually the 
case. Reading now off the weight of each state (sf | G Sm, and multiplying by 
we deduce the value of the matrix element (sf |T(g)|si). 

A few practical remarks are in order. Since all weights are polynomials in q the 
computations are done using symbolic algebra, representing the integer coefficients of 
each power of q in an array. These integers may in principle be very large, and thus 
cause overfiow in a standard 32 bit integer arithmetic. This problem can be easily 
remedied by using modular arithmetic, i.e. computing all coefficients modulo differ- 
ent primes, and retrieving the complete result from the Chinese remainder theorem. 
However, for the system sizes used in this work it turned out that 32 bits are actually 
sufficient. 

As usually with transfer matrix methods, the main obstacle for going to large 
system sizes is due to memory limitations. Our algorithm for m = 13 used some 500 
MBytes of memory, and although one might have gone to slightly larger systems on 
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a bigger computer, we refrained from doing so because of the restrictions imposed 
by the subsequent analysis (see Section H). The efficiency of our algorithms can be 
appreciated by noting that the calculations for m = 8, the largest system size treated 
pT| , were done in less than ten seconds. In fact we have reproduced all results 



m 



given in Ref. [O] for the transfer matrices and the boundary condition vectors (see 



below) with m < 8; besides validating Ref. |11|, this is of course also a nontrivial 
check on the correctness of our programs. 

Apart from the transfer matrix we are also interested in computing the zero- 
temperature partition functions (i.e., the chromatic polynomials) P{mp x rip) for 
strips of a finite length n and with free boundary conditions in the longitudinal 



direction. As explained in Section 6] this can be achieved by sandwiching powers 



of T{q) between suitable initial and final vectors, 

P{mpXnp) = {vi\T{qr-'\v,). (2.4) 

The initial vector \v]) is simply the unit vector defined by assigning a weight one 
to the unique state in which each of the m points belongs to a distinct cluster, and 
zero to all other states. The final vector (ff | corresponds to adding the last row of the 
lattice, followed by a trace over the final states (sf | in which each state is weighted by 
qki{si\)_ jg found by a procedure very similar to the one used to compute the matrix 
elements (sf |T(g)|si). The only difference is that, in the last step, we take the trace 
over the (sf|, rather than projecting on each of them and multiplying by g~'^(^^fl). 



3 Numerical Results for the Square-Lattice Chro- 
matic Polynomial: 

Periodic Transverse Boundary Conditions 

We have computed the transfer matrix T{q) and the limiting curves B for square- 
lattice strips of widths 9 < L^. < 13 with periodic boundary conditions in the trans- 
verse direction. One direct check of our results is provided by the (trivial) identity 



Z{mp X Up) = Z{np x mp) 



(3.1) 



The chromatic polynomials Z{np x mp) with n = 2,3 are well-known p2| , [29[] . We 



have used both the resultant and the direct-search methods to obtain the limiting 
curves (see |]11| for details). To compute the isolated limiting points we define the 
matrix D[q): 

Pnx2 
Pnx3 



( P 



D{q) 



nxl 
Pnx2 



PnxM 
P 

J. ri 



\ 



nxM+l 



\ PnxM ^nxM+l 



Pr, 



(3.2) 



<2M-1 j 
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where M is the dimension of the corresponding transfer matrix. Then it can be proved 
that 

M 

detD(g) = J]afc J] (A, - A,)^. (3.3) 

k=l l<i<j<M 



Thus, all the zeros of the amplitudes ak are also zeros of det D{q) (see also [|TT], Section 
2.2]). 



3.1 Lx — 9p 

The transfer matrix is 22-dimensional; it can be found in the Mathematica file 
transf er2 .m that is available with the electronic version of this paper in the cond-mat 
archive at [tittp : / /www . lanl . gov| . We have checked that none of the amplitudes 
vanishes identically. 

We have computed the limiting curve B using the direct-search method. In addi- 
tion, we were able to compute the resultant at 6* = (but unfortunately not for all 
6). Some of its zeros are precisely the endpoints of the limiting curve, which is why 
our results for the endpoints are more accurate in this section than in the following 
ones.Q The limiting curve contains nine disconnected pieces (see Figure [1|). One of 
them crosses the real axis at g ~ 2.8505722621. There are 20 endpoints: 



-0.44699181 ± 1.24252295 i 
1.14912991 ±2.61264326 i 
2.36507014 ± 2.08719033 i 
2.83939211 ± 1.19488272 i 
2.84682005 ± 0.40527862 i 
3.01123194 ± 1.04259193 i 



1.15312070 ± 2.61981372 i 
2.38763742 ± 2.08307229 z 
2.83111413 ± 1.16708976 i 
2.85144522 ± 0.40515929 i 



(3.4a) 
(3.4b) 
(3.4c) 
(3.4d) 
(3.4e) 
(3.4f) 



There are eight small gaps defined by the endpoints given in Eqs. (3.4b ^3^) . Finally, 
there are two T-points at g ~ 2.854 ± 1.123 i. 

We have also computed the polynomial det D{q) for this strip. It can be written 



as 



det D{q) = q'\q-iy\q-2y\q' -3q + iy%q~3y\q' ~5q' + Qq~iy 

(g2 - 4g + 2)^(g3 - 6?^ + 9q - l){q'^ - 5q + h)P{qf, (3.5) 

where P{q) is a polynomial of degree 1170 with integer coefficients (and not fac- 
torizable over the integers). The full expression of (|3.5|) is included in the Math- 
ematica file transf er2 .m. The pre-factors appearing in ( p.5| ) correspond to the 
first nine minimal polynomials Pri(q') given in [0, Table 1]. In particular, this 



^ The zeros of the resultant at = correspond to values of q such that there are two (or more) 
equal eigenvalues. However, we are only interested in the values of q associated to dominant equal 
eigenvalues. 
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means that det D{q) = for the first nine Beraha numbers B2, . . ■ , Biq, in agree- 
ment with Conjecture 7.1 of |Tl|. Indeed, the first three Beraha numbers 0,1,2 



are trivial zeros as all the amplitudes vanish at those values. Thus, they are iso- 
lated limiting points. The non-trivial isolated limiting points are the Beraha number 
B5 = (3 + y5)/2, and the complex conjugate pairs q ^ 1.3290624037±2.5498775371 i 
and q ^ 2.5285781385 ± 1.5750286926 i (see Figure [l|). 

The convergence to the non-trivial real zero is quite fast (actually, exponential) 
as shown in Table |l|. The largest real zero in Table |l| converges at an approximate 
1/n rate to the value go ~ 2.8505722621. 

We have also checked Conjectures 7.2 and 7.3 of |jll|, namely that det D{q) > 



for all q = Bk with k > 10. In particular, we have checked this result up to A; = 50. 
We expect that the same result holds also for > 50 as f{k) = det D{q = B^) is 
clearly a increasing monotonic function of k. 

3.2 = lOp 

The transfer matrix is 51-dimensional; it can be found in the Mathematica file 
transfer2.m. We have checked that none of the amplitudes vanishes identically. 
Unfortunately we were unable to compute the resultant even at 6* = 0, so we used a 
pure direct-search method. 

The limiting curve B contains ten connected pieces. None of them crosses the 
real axis. The points closest to that axis are q ~ 2.82813 ± 0.00097 2. We find the 
following 22 approximate endpoints: 

q ^ -0.48079 ± 1.13785 i (3.6a) 

q ^ 0.94694 ± 2.61250 i; q^ 0.94561 ± 2.60984 i (3.6b) 

q ^ 2.18917 ± 2.26882 i; q^ 2.19871 ± 2.26845 z (3.6c) 

q ^ 2.76451 ± 1.49846 i; q^ 2.77167 ± 1.48200 i (3.6d) 

q ^ 2.82813 ±0.00097i (3.6e) 

q ^ 2.86581 ±0.87778 q^ 2.87243 ± 0.87284 z (3.6f) 

q ^ 3.21691 ±0.52738 i (3.6g) 



There are nine small gaps defined by the endpoints (3.6b- 3]6|). Finally, there is a 
T-point close to g ^ 2.874 ± 0.941 i. 

In this case we have been unable to obtain the polynomial det D{q) due to lack 
of computer memory. By inspection of Figure ^ we expect four real limiting points 
at the first four Behara numbers q = 0,1,2, B^. The first two limiting points (i.e., 
q = 0, 1) are trivial ones, as all the amplitudes vanish. We have computed the 
amplitudes corresponding to g = 2, B^ and, in both cases, the dominant amplitude 
a* vanishes. The convergence to the non-trivial real zeros g = 2, i?5 is quite fast as 
shown in Table 0. There are also three pairs of complex conjugate isolated limiting 
points: g ^ 2.7535938496 ± 0.9533568270 i, g ^ 2.3640435205 ± 1.8567782883 and 
g ^ 1.0488718011 ±2.5727097899 i. We have computed the amplitudes corresponding 
to these values of g. The dominant amplitude is very small for the first two pairs (the 
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absolute value of that amplitude is respectively ~ 4.09 x 10"'^^, and ~ 6.4 x 10~'^^). 
The dominant amplitude for the last pair is not that small ~ 3.03 x 10~^). 

We beheve there is an isolated limiting point nearby, but we cannot guess a better 
estimate as the convergence is still dominated by the 1/n ratio coming from the 
regular points of B. This phenomenon also occurs for = 9p (see the upper isolated 
limiting point in Figure p. 

We have also computed the amplitudes ak at the first Beraha numbers Bn with 
n < 22. For B2, . . . , i?5 the dominant amplitude vanishes as explained above. For 
Bq, . . . , Bii we find at least one vanishing amplitude, but none of them is the dominant 
one. Thus, none of these Beraha numbers is an isolated limiting point. Finally, for 
B12, . . . , B22 all the amplitudes are nonzero. Furthermore, the product of all the 
amplitudes for these Beraha numbers is positive; this implies that det D{q) > for 
q = Bi2,...,B22 because of Eq. (p.3|) . These results support Conjectures 7.1, 7.2, 



and 7.3 of ITT 



3.3 Lx = lip 

The transfer matrix is 9 5- dimensional; it can be found in the Mathematica file 
transfer2.m. We have checked that none of the amplitudes vanishes identically. 

The limiting curve crosses the real axis at g ~ 2.8900930977. We find the following 
20 approximate endpoints: 

q ^ -0.50104 ± 1.05084 i (3.7a) 



q ^ 0.77281 ±2.58699 i 

q ^ 2.01541 ± 2.39925 i 

q ^ 2.66452 ± 1.73992 i 

q ^ 2.88992 ± 0.35265 i 



q^ 0.77329 ± 2.58795 z (3.7b) 

q^ 2.01168 ± 2.39892 i (3.7c) 

q^ 2.67029 ± 1.73369 i (3.7d) 

q^ 2.89068 ± 0.35252 z (3.7e) 



q ^ 3.09853 ± 0.85628 i (3.7f) 

The endpoints (3.7b- pl7^ define eight small gaps. Finally, there is a pair of complex 
conjugate T-points close to g ~ 2.902 ± 1.015 i 

By inspection of Figure ^ we expect four real limiting points at the first four 
Beraha numbers q = 0,1, 2, £5. The first three values (i.e., q = 0,1,2) corre- 
spond to trivial isolated limiting points as all the amplitudes vanish. The fourth 
value q = B^ has a vanishing dominant amplitude; thus, it is a true isolated lim- 
iting point. The convergence to g = -B5 is quite fast as shown in Table |1[ The 
largest real zero in Table |l] converges at an approximate 1/n rate to the value 
go ~ 2.8900930977. There are also three pairs of complex conjugate isolated lim- 
iting points: g ^ 2.6642077467 ± 1.2375206436 i, q ^ 2.1718870063 ± 2.0937620001 i, 
and g ~ 0.8331263199 ± 2.5614743076 2. The value of the dominant amplitude van- 
ishes (within our numerical precision) for the first two pairs (namely, ~ 10~^°). 
However, the last pair is close, but not quite equal, to the real isolated limiting point. 
In this case we have ^ 4.70 x 10~^. We believe that there is an isolated limiting 
point nearby, but that the convergence to this point is still dominated by the 1/n 
rate coming from the nearby regular points in B. 



10 



We have also checked that Conjectures 7.1, 7.2, and 7.3 of [11] hold for this strip: 
there is at least one vanishing amplitude for q = B2, . . . , -B12, while all the amplitudes 
are nonzero for q = B13, . . . , B22- Only ioi q = B2, . . . , the dominant amplitude 
vanishes; thus, these are the only real isolated limiting points. Finally, the product 
Ilfc "fc > for Q = Bi3, . . . , -B22, implying that det D{q) > 0. 



3.4 = 12 



The transfer matrix is 232-dimensional; it can be found in the Mathematica file 
transfer2.m. We have checked that none of the amplitudes vanish identically. We 
have not been able to compute the limiting curve B for this the transfer matrix 

is too large. However, we could obtain the chromatic polynomials for square-lattice 
strips with aspect ratios between 1 and 8 (i.e., Pi2xi2, -Pi2x24, • • • , -Pi2x96)- As shown 
in Figure ^, the zeros converge to a limiting curve very similar to the ones obtained 
in the previous cases. Via the direct-search method (applied to the limiting curve) 
we have been able to obtain an estimate of the point of the limiting curve B which is 
closest to the real axis, reading q ^ 2.874356 ± 0.00037 z. 

By inspection of Figure ^ we conclude that there are four real isolated limiting 
points at the first four Beraha numbers. The first two of those points (i.e., g = 0, 1) 
are trivial as they correspond to the vanishing of all the amplitudes. The convergence 
to the other two (non-trivial) real limiting points is quite fast as can be observed 
from Table [^. We also find that the dominant amplitudes for q = 2,B^ are very 
small. This evidence supports our belief that g = 2 and B^ are true isolated limiting 
points. We also find by inspection three pairs of complex conjugate isolated limiting 
points: q ^ 1.9629324150 ± 2.2751179548 i, q ^ 2.5667183607 ± 1.4761470034 and 
q ~ 2.7956875974 ± 0.7913632912 z. The values of \ai,\ are very small for these values 
of q: 1.11 X 10-l^ 5.29 x lO'^o, and 1.22 x 10"^*^ respectively. 



We have also checked that Conjectures 7.1, 7.2, and 7.3 of [jlj] hold for this 
strip: there is at least one vanishing amplitude for q = B2, . . . , -B13, whereas all the 
amplitudes are nonzero for q = i?i4, . . . , i?22- Furthermore, for these latter values 
Hfcafc > 0, thus det D{q) > ( ^31) . 



3.5 Lx — 13p 

The transfer matrix is 498-dimensional; it can be found in the Mathematica file 
transfer2.m. We have checked that none of the amplitudes vanish identically. We 
have not been able to compute the limiting curve B for this case as the transfer matrix 
is too large. However, we could obtain the chromatic polynomials for square-lattice 
strips with aspect ratios between 1 and 5 (i.e., -Pisxis, -Pisxes)- As shown in 
Figure ^, the zeros converge to a limiting curve very similar to the ones obtained in 
the previous cases. Even though the computation of the full limiting curve is too time 
consuming, it is very interesting to have an estimate for the value where B crosses 
the real axis. Our result is go ~ 2.9161885031, using the direct search method. 

By inspection of Figure ^ we conjecture that there are 4 real isolated limiting 
points corresponding to the first four Beraha numbers g = 0, 1, 2, i?5. The first three 
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are trivial zeros as all the amplitudes vanish. From Table ^ we observe that the 
convergence to g = -B5 is very fast. We have also computed the dominant am- 
plitude for q = B^. This indeed vanishes, thus g = -B5 is a true isolated limit- 
ing point. The largest real zero in Table |l] converges at an approximate 1/n rate 
to the value qo ~ 2.9161885031. We also find empirically three pairs of complex 
conjugate zeros that are likely to be isolated limiting points. As the convergence 
to those points is exponentially fast, the estimates obtained from the chromatic 
zeros for the lattice 13p x 53f are expected to be close enough to the true val- 
ues: q ^ 1.7591342818 ± 2.3988734769 i, q ^ 2.4614147948 ± 1.6813861639 i, and 
q ^ 2.7323074066 ± 1.0328078058 Indeed, we obtain very small dominant ampli- 
tudes for all three values of g, namely I I ^ 7.48x10"^^, 1.30 x 10~^^, and 1.18 x 10~^^ 
respectively. 



We have also checked that Conjectures 7.1, 7.2, and 7.3 of [jl^] hold for this 
strip: there is at least one vanishing amplitude for g = i?2, • • • , -B14, whereas all the 
amplitudes are nonzero for g = B15, . . . , B22- Furthermore, for these latter values 
Hfc ctk > 0, thus det D(g) > 



4 Conclusions 

We have computed the transfer matrix of the zero-temperature Potts antiferro- 
magnet defined on square-lattice strips of widths 9 < m < 13 with cylindrical bound- 
ary conditions; the corresponding result with free boundary conditions are available 
from the authors upon request. For m = 9, 10, 11 we have been able to compute the 
limiting curves of the partition function zeros using the direct-search method. All 
of them have a similar qualitative shape (see Figure ^) although we observe slight 
differences between the limiting curves with even width and those with odd width. 

In Table |^ we summarize the basic features of the limiting curves. It is worth 
noting that the number of connected components and endpoints increases with the 
strip width m. However, the gaps between those endpoints narrow as m increases. 
Thus, the thermodynamic limit m — 00 seems to be rather complicated in this model. 
This may be the cause the solution to this model is still lacking. On the other hand, 
the thermodynamic limit for the triangular-lattice case is smoother |3^ and the model 
is solvable . 



The differences between even- and odd-width strips is manifest when we try to 
define go. Only the curves corresponding to odd width cross the real axis (see Table ||). 
Thus, the value of go can be defined for those strips. We also observe that the value 
of go is monotonically increasing with the width m. Thus we conjecture that 

Conjecture 4.1 The limiting curves for square-lattice strips with cylindrical bound- 
ary conditions and odd width m cross the real q-axis at a point qQ{m). The function 
go(m) is a monotonically increasing function ofm. 



Our best estimate for go (assuming that Conjecture |]T] holds) is 

go(sq) > 2.9161885031 . (4.1) 
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This estimate comes from the strip with = 13p. This value of go is surprisingly 
close to the expected value for gc(sq) = 3 for this lattice. (It is closer to qc than in case 
of the triangular lattice where qQ{iii) ~ 3.81967 and gc(tri) = 4.) It might happen 
that for the square lattice go = Qc- If this were the case, the limiting curve B when 
m ^ oo would be qualitatively different from the triangular-lattice limiting curve 
computed by Baxter |^|. In particular, instead of the three domains of definition of 



the free energy on the positive real g-axis for the triangular-lattice case, we would 
have only two such domains. 

In Table @ we also observe that for even strip widths m the limiting curve does 
not cross the real axis. We can nevertheless define go to be the closest endpoint to the 
real g-axis (with positive imaginary part). The real part of this quantity seems to be 
also a monotonically increasing function of the lattice width, while the imaginary part 
goes rather quickly to zero. Assuming the monotonicity of Re go our best estimate 
for go would be go(sq) ^ 2.87436 which is smaller than ( [4.11 ). Let us now consider the 
quantity Imgo and try to fit it to a power-law Ansatz. We obtain Imgo ~ 26.4m~^'^^, 
indicating that the true behavior might in fact be exponential. Thus, we tried an 
exponential Ansatz Imgo = ^ x . We obtain a very good fit including all the data: 
A = 0.1177(39) and B = 0.6188(26) with = 0.019 (1 degree of freedom (DF), 
confidence level (level) = 89%).^ The similarity between B^^ = 1.6160(67) and the 
golden ratio r = bI^^ = {1 + ^)/2^ 1.6180339887 is striking. One may speculate 
in a possible connection between this appearance of B^ and its role as the largest real 
isolated limiting point in the m —>■ oo limit. Thus, we state the conjecture 



Conjecture 4.2 The limiting curves for square-lattice strips with cylindrical bound- 
ary conditions and even width m do not cross the real q-axis. Let qo{m) be the closest 
point of the limiting curve to the real q-axis. Then, Rego(m) is a monotonically in- 
creasing function of m, and Imgo(m) ~ B^ ^^"^ = where r = (1 + V5)/2 is the 
golden ratio. 

Let us go back to the quantity go for lattice strips with odd m. We can try to fit 
the data to the power-law Ansatz A + Bm~^. In all the fits we impose a lower cutoff 
?7T. > mmin to the data and we study the effect of the cutoff on the estimates A, B 
and A and on the x^. As we increase mmin the estimate A decreases from 3.12688(2) 
down to 3.02799(5); the estimate -B increases from 1.49601(4) to 2.8011(19); and 
the exponent A increases from 0.76651(40) to 1.25578(42). The is very poor: 
8.69 X 10^ and 1.75 x 10^ for m^^^ = 5 and 7 respectively. It seems that the estimate 
for A tends to 3 from above. As we know that the true value for A should be bounded 
from above by qc = 3, we can try to fit the data to the Ansatz: 3 — go = Bm~^. 
Again, we do not find any stable fit: as mmin increases from 5 to 11, the estimate 
B monotonicly goes from 2.47978(4) to 5.3807(12), and the exponent A grows from 
1.28668(1) to 1.62265(9). Again the is poor: it goes from 1.72 x 10*^ for mmin = 5 
to 5.22 X 10^ for m^i^ = 9. 



^ We have fitted the data shown in Table g; the error bars are assigned to be one unit in the least 
significative digit. 
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We can play the same game with the quantity Re go for the lattice strips of even 
width m. The fits to the Ansatz A + Bm~^ are similar to the former case: as m 
is increased, A decreases from 3.22832(3) down to 3.02285(28); —B increases from 
2.95545(6) to 5.970(19); and A increases from 0.87657(4) to 1.4866(20). The is 
poor: 5.19 x 10^ and 5.30 x 10'^ for mmin = 4 and 6 respectively. Again the estimate 
for A seems to tend to 3 from above, so we can try the next Ansatz 3 — Re go = Bm~^. 
These fits are not very stable either: as mmin increases from 6 to 10, -B monotonicly 
goes from 6.11369(19) to 8.989(12), and A grows from 1.54057(2) to 1.71849(54). The 

is also poor: it goes from 7.02 x 10^ for mmin = 6 to 8.45 x 10'^ for m^in = 8. 

Our data thus suggest that, within our numerical precision, go(sq) = gc(sq) = 3. 
We can state this finding conjecture 

Conjecture 4.3 For the zero-temperature square-lattice Potts antiferromagnet with 
cylindrical boundary conditions, Qq = qc = 3. 



The various limiting curves B are most similar for intermediate values of Reg 
(namely, 0.6 < Reg < 2.4). The finite-size effects are more apparent for g close to 
g = and g = qc- This is expected from what Baxter |]^ found for the triangular 



lattice: the density of points around g = and g = gc was smaller than for the other 
regions in the complex g-plane; thus, the convergence is expected to be slower around 
those two points. 

In particular, we expect that eventually the rightmost branches of each limiting 
curve (see Figure ||) will curve back and decrease towards the expected value gc = 3. 
Indeed, this phenomenon of "overshooting" seems to be characteristic of cylindric 
boundary conditions; comparison with the limiting curves B found in Ref. |TT| for 



free boundary conditions shows that in this case the approach towards the limiting 
value gc = 3 is monotonic. A similar remark holds true for the Potts antiferromagnet 
on the triangular lattice ||33[. On the other hand, the leftmost endpoints of the 
limiting curves B are expected to curve back and converge to the point g = as in 
the triangular-lattice case. 



We have also checked Conjectures 7.1, 7.2, and 7.3 of Ref. |jTT|, Section 7]. In 
particular, they state that for a square-lattice strip of width L with free or cylindrical 
boundary conditions, 

1. At each Beraha number q = B2, . . . , -Bl+i there is at least one vanishing ampli- 
tude afc(g); 

2. For all q = B^ with k > L + 1 none of the amplitudes afc(g) vanishes; and 

3. det D{q) > for all q = Bk with k> L + 1. 

We have shown that Conjectures 7.2 and 7.3 hold up to /c = 50 for L = 9p, and up to 
k = 22 for lOp < L < 13p. In all cases, we have found empirically that the number of 
vanishing amplitudes decreases from q = B2 to q = i?L+i- In this latter case there is 
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only one vanishing amplitude and its corresponding eigenvalue is given by A = (— 1)^. 
We have checked that this is so for all L < 13p and < 8f.[] 
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Lattice 



9p X 9f 
9p X 18p 
9p X 27p 
9p X 36p 
9p X 45p 
9p X 54p 
9p X 63f 
9p X 72f 
9p X 81f 
9p X 90f 



lOp X 
lOp X 
lOp X 

lOp X 

lOp X 
lOp X 
lOp X 
lOp X 
lOp X 
lOp X 



IOf 
20f 
30f 
40f 

50f 
60f 
70f 
80f 
90f 
lOOp 



lip X 
lip X 
lip X 
lip X 
lip X 
lip X 
lip X 
lip X 
lip X 
lip X 



lip 

22f 
33f 
44f 
55f 
66f 
77p 
88f 
99f 
IIOf 



12p X 12p 

24f 
36f 
48f 
60f 
72f 
84f 
96f 



12p 
12p 
12p 
12p 
12p 
12p 
12p 



13p X 
13p X 
13p X 
13p X 
13p X 



13f 
26f 
39f 
52f 
65p 



Beraha 



3rd Zero 



2.000000000000 
2.000000000000 
2.000000000000 
2.000000000000 

2.000000000000 
2.000000000000 
2.000000000000 
2.000000000000 
2.000000000000 
2.000000000000 



2.000000000000 

2.000000000000 
2.000000000000 
2.000000000000 
2.000000000000 
2.000000000000 
2.000000000000 
2.000000000000 



4th Zero 



2.620608391171 
2.618033988527 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 



2.614232442467 
2.618033987905 
2.618033988750 
2.618033988750 

2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 



2.618035720465 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 



2.618032077436 

2.6180.33988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 



2.618033988926 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 



2.618033988750 



5th Zero 



2.660519126718 
2.791347005408 
2.814819876147 
2.824928798467 
2.830572841537 



2.731088786113 
2.836691210598 
2.857656644942 
2.866767905070 
2.871875414535 



2.775231591227 
2.867878041004 
2.886773116380 



Table 1 : Real zeros of the chromatic polynomials of finite square- lattice strips with 
periodic boundary conditions in the transverse direction and free boundary conditions 
in the longitudinal direction, to 12 decimal places. A blank means that the zero in 
question is absent. The first two real zeros g = 0, 1 are exact on all lattices; the third 
real zero g = 2 is exact on all lattices of odd width. "Beraha" indicates the Beraha 
numbers B4 = 2 and = (3 + \/5)/2. 
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Eigenvalue-Crossing Curves B 


Isolated Points 


Lattice 


#C 


#E 


#T 


#D 


#ER 


min Re q 


qo 


max Re q 


#RI 


#CI 


3p 























3 





4p 


3 


8 





1 





0.709803 


[2.253370,2.351688] 


2.995331 


3 





5p 


5 


10 











0.165021 


2.691684 


2.691684 


4 





6p 


3 


10 


2 


1 





-0.131889 


[2.608943, 2.613228] 


3.171192 


3 


1 


7p 


7 


16 


2 








-0.296250 


2.788378 


2.839016 


4 


1 


8p 


6 


16 


4 








-0.390864 


2.751531 ±0.002531 r 


3.211132 


4 


2 


9p 


9 


20 


2 








-0.446992 


2.850572 


3.011232 


4 


4 


lOp 


lot 


22t 


2t 


ot 


Qt 


-0.48079 


2.82813 ±0.00097 i*. 


3.21691 


4 


6t 


lip 


9t 


20t 


2t 


ot 


Qt 


-0.50104 


2.890093 


3.09853 


4 


6t 


12p 














2.87436 ± 0.00037 i* 




4 


6t 


13p 














2.916189 




4 


6t 



Table 2: Summary of qualitative results for the eigenvalue-crossing curves B and 
for the isolated limiting points of zeros. For each square-lattice strip considered in 
this paper, we give the number of connected components of B C), the number 
of endpoints E), the number of T points T), the number of double points 
D), and the number of enclosed regions ER); we give the minimum value of Reg 
on B, the value(s) go where B intersects the real axis (* denotes an almost-crossing), 
and the maximum value of Reg on B. We also report the number of real isolated 
limiting points of zeros RI) [which are always successive Beraha numbers B2, -B3, 
. . .] and the number of complex conjugate pairs of isolated limiting points CI). 
The symbol ^ indicates uncertain results. The results for L < 8p are included for 



comparison; they are taken from Ref. [11 
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Zeros sq lattice L = 9 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




I I I I I I I I I I I I I I I I I I I 

12 3 

Re(q) 

Figure 1: Zeros of the partition function of the g-state Potts antiferromagnet on 
the square lattices 9p x 45i? (squares), 9p x 90f (circles) and 9p x oop (solid line). 
The isolated limiting zeros are depicted by a x. The limiting curve was computed 
using the direct-search method, except the endpoints that were computed using the 
resultant method. 
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Zeros sq lattice = 10 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




J I I I I I I I I I I I I I I I I I L 



12 3 

Re(q) 

Figure 2: Zeros of the partition function of the g-state Potts antiferro magnet on the 
square lattices lOp x 50^ (squares), lOp x lOOp (circles) and lOp x oop (solid line). 
The isolated limiting zeros are depicted by a x. The limiting curve was computed 
using the direct-search method. 
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Zeros sq lattice L = 



lip 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




J I I I I I I I I I I I I I I I I I L 



12 3 

Re(q) 

Figure 3: Zeros of the partition function of the g-state Potts antiferro magnet on the 
square lattices lip x 55^ (squares), lip x llOp (circles) and lip x oop (solid line). 
The isolated limiting zeros are depicted by a x. The limiting curve was computed 
using the direct-search method. 
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Zeros sq lattice = 12 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




J I I I I I I I I I I I I I I I I I L 



12 3 

Re(q) 

Figure 4: Zeros of the partition function of the g-state Potts antiferromagnet on the 
square lattices 12p x 60 p (squares) and 12p x 96p (circles). The isolated limiting 
zeros are depicted by a x. 
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Zeros sq lattice = 13 
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Figure 5: Zeros of the partition function of the g-state Potts antiferromagnet on the 
square lattices 13p x 52 p (squares) and 13p x 65p (circles). The isolated limiting 
zeros are depicted by a x. 
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Limiting Curves Square Lattice 
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Figure 6: Limiting curves for the square-lattice strips Lp x oop with 8 < L < 11. 
The curve for L = 8 was obtained in 111 111. 
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